Model of bound interface dynamics for coupled magnetic domain walls 
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A domain wall in a ferromagnetic system will move under the action of an external magnetic 
field. Ultrathin Co layers sandwiched between Pt have been shown to be a suitable experimental 
realization of a weakly disordered 2D medium in which to study the dynamics of ID interfaces 
(magnetic domain walls). The behavior of these systems is encapsulated in the velocity-field response 
v(H) of the domain walls. In a recent paper [P.J. Metaxas et ai, Phys. Rev. Lett. 104, 237206 
(2010)] we studied the effect of ferromagnetic coupling between two such ultrathin layers, each 
exhibiting different v(H) characteristics. The main result was the existence of bound states over 
finite-width field ranges, wherein walls in the two layers moved together at the same speed. Here, 
we discuss in detail the theory of domain wall dynamics in coupled systems. In particular, we show 
that a bound creep state is expected for vanishing H and we give the analytical, parameter free 
expression for its velocity which agrees well with experimental results. 

PACS numbers: 75.78.Fg, 75.60.Ch, 75.70.Cn 



I. INTRODUCTION 

A number of physical phenomena involve elastic inter- 
faces moving through disordered media. These phenom- 
ena range from domain wall motion in ferromagnets^— , 
ferroelectrics 5 and multifcrroics 6 to wetting 7 as well as 
vortex motion in high-T c superconductors^. The theo- 
retical frameworks^— developed to model elastic inter- 
face dynamics are therefore highly relevant for a num- 
ber of real world processes which are of interest both for 
their fundamental properties and eventual applications. 
Indeed, theoretical studies of single interface dynamics 
and statics have revealed a lot of interesting physics with 
predictions of universality, in-depth studies of dynamic 
and static critical exponents 9 - - — and the development of 
now well-known interface growth equations^. Magnetic 
systems in particular have been an ideal testing ground 
for these theories^— ^ since these systems can be easily 
probed and manipulated. 

A relatively recent theoretical, and more recently, ex- 
perimental, playground has been developing concern- 
ing the physics of interacting interfaces in 2D sys- 
tems. Theoretically, this problem has been studied 
via modified growth equation s 16 ' 17 , Monte Carlo mod- 
eling of repulsive or non-interacting interface o 18 ' 19 and 
scaling arguments 2 ^. Quasi-2D experimental realiza- 
tions of systems containing coupled interfaces have also 
been conceived, ranging from interacting fluid fronts^ 1 - to 
repulsive— and attractiv e) 22 ' 31 magnetic domain walls. 

Our work on field-driven, attractively coupled domain 
walls 22 has been carried out on a system consisting of 
two physically separate, but magnetically coupled^ 3 -, ul- 



a) 



Pt (3.5 nm) 



Co (0.8 nm) 



Pt (3 nm) 



Co (0.5 nm) 



Pt (4.5 nm) 



Domain wall 




J>0 



Domain wall 



b) 




c) 




Coupled domain walls 



Single domain wall 



FIG. 1. (a) Using magnetic layers coupled via an interlayer 
interaction energy, J, to create a model system for studying 
bound domain walls. Measurements of domain wall motion 
in this system consisted either of (b) coupled domain walls 
(boundaries of domains which are statically aligned in zero 
field) or (c) single domain walls (boundaries of domains ex- 
isting in the hard layer only). 



trathin ferromagnetic Co layers [Fig. [Ha)]. The ferro- 
magnetic coupling tends to align the magnetization in 
the two layers. Therefore, if a domain wall is present 
in each Co layer [eg. Fig. QJb)], the ferromagnetic cou- 
pling will tend to align them, acting as an attractive in- 
teraction between the walls. This attraction not only 
favors a static domain wall alignment in zero-field, but 
can also stabilize the aligned state dynamically under 
an applied fielc-2 2 -. In this case, walls in the two layers 
are dynamically bound and move together at a common, 
unique velocity, despite each wall having different intrin- 
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sic velocity-field responses. These differing velocity-field 
responses however do mean that dynamic domain wall 
binding can occur only over field ranges in which wall ve- 
locities in each layer are sufficiently close, placing a limit 
on the fields for which bound motion can occur. Un- 
til this work, studies of pairs of interacting interfaces in 
quasi-2D systems had been mostly carried out in single 
media. While it was already thought that domain walls 
in strongly coupled layers moved together— ~—, this was 
the first study wherein both dynamically bound domain 
walls and transitions between bound and unbound dy- 
namics were directly evidenced. 

In this article, we discuss in detail a theoretical descrip- 
tion of bound domain wall motion. The paper is outlined 
as follows. In Sec II we briefly give some details about 
the model system. In Sec. Ill we analyze how domain 
wall speed is affected by interlayer coupling and in Sec. 
IV we study analytically the bound state regimes and 
discuss the agreement between theory and experiment. 
A short conclusion follows. 



II. COUPLED ULTRATHIN MAGNETIC 
LAYERS 



The experimental system shown in Fig. HJa) is a mag- 
netic multilayer consisting of two ultrathin Co layers: a 
magnetically hard 0.8 nm layer (layer 1) and a softer 
0.5 nm layer (layer 2). The layers are ferromagneti- 
cally coupleo-2, (coupling energy J > 0) across a 3 nm 
thick Pt spacer. Seed and capping Pt layers ensure an 
out-of-plane magnetic anisotropy within the Co layers. 
Pt/Co-based films are now considered good experimen- 
tal realizations of a weakly disordered, ferromagnetic 
2D Ising system, due to their anisotropy-induced out-of- 
plane magnetization, narrow domain walls and intrinsic 
structural disorder^ 3 -. This disorder has a major role in 
determining the velocity response v(H ) of a domain wall 
to an external field H, applied perpendicular to the film 
plane. 

Two types of domain wall velocity measurements were 
carried out based upon the two domain (wall) types 
which could be nucleated within the multilayer. Both 
types of wall could be propagated under field to deter- 
mine their velocity- field responses using a quasi-static 
magneto-optical method 3 ^-. (1) Coupled domain walls 
are the boundaries of domains existing in both layers 
which, in zero field, are aligned spatially with their mag- 
netizations pointing in the same direction, as shown in 
Fig- Hlb). Under field, and depending on the field ampli- 
tude, they can move together, in a dynamically bound 
state, or separately. (2) Single domain walls are the 
boundaries of domains existing in the hard layer only, 
as illustrated in Fig. [He). Measurements of these do- 
main walls yield a reference velocity and a determination 
of the interlayer coupling. 



III. FROM ISOLATED TO COUPLED AND 
BOUND DOMAIN WALL DYNAMICS 

Here we analyze field-velocity responses of: i) a single 
domain wall in an isolated magnetic layer, ii) a single do- 
main wall in one magnetic layer coupled to a second, sat- 
urated magnetic layer, and iii) two coupled domain walls, 
one in the hard layer and the other in the soft layer. Do- 
main walls will be approximated as straight lines, whose 
position is given by a single number. We begin with a 
single wall located at x — x w in an isolated ultrathin Co 
layer [see Fig. [2fa)]. The Co layer is positively magne- 
tized for x < x w and negatively magnetized for x > x w . 
The application of an external field H > drives the 
wall to the right, with the wall acquiring a positive ve- 
locity v(H) — dx w /dt. Experimental results 3 obtained 
for domain wall motion in Pt/Co(0.5-0.8 nm)/Pt films 
show that v(H) is characterized by two distinct regimes 
at room temperature (creep and flow) which were theo- 
retically predicted^ and are sketched in the schematic of 
Fig.[5{b). Domain walls exhibit flow motion at high fields 
for which v oc H. However, below a layer- dependent crit- 
ical depinning field H dBp (generally on the order of a few 
hundred Oersted 3 -), disorder- induced pinning effects be- 
come significant and the walls exhibit thermally activated 
creepi. Within this latter regime, v(H) has the form 



v(H) — vq exp 



U c [H da 



k R T 



H 



1/4' 



(1) 



where the exponential factor Uc/ksT is the ratio be- 
tween the typical pinning energy and the thermal energy. 
The exponent 1/4 is a universal exponent, characteristic 
of the dynamics of a one dimensional interface in a 2D 
weakly disordered medium. 

Films with different thicknesses have different micro- 
scopic parameters and disorder strengths. As a result, 
they have different v(H) characteristics^, as attested by 
the experimental velocity-field curves for domain walls in 
the two layers in the absence of coupling [Fig.[2][c)]. How- 
ever, pairs of such curves often intersect at two points: 
H = and H = H* (H* ~ 860 Oe for our system). 
The first crossing point is universal, because the velocity 
v(f) of any isolated interface in response to a generalized 
force / (here, / = H) is always expected to vanish for 
vanishing /. The second crossing point is less trivial and 
arises because domain walls in thicker Co layers gener- 
ally have a lower creep velocity but a higher flow velocity 
than walls in thinner layers. 

In the remainder of this article, we shall use v\(H) 
to refer to the domain wall velocity in the hard layer 
and V2(H) to that in the soft layer. If the two films are 
not coupled, it is clear that walls will propagate indepen- 
dently, with v 2 > v\ for H < H* and V\ > v 2 for H > H* 
[Fig. HJc)]. The question we are now going to consider is 
the following: What is the effect of interlayer coupling on 
domain wall velocities v\^{H) and domain wall binding 
phenomena? 
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FIG. 2. (a) The average domain wall position is denoted x w . 
(b) At finite temperature, walls exhibit a low field, thermally 
activated creep regime and a high field, dissipation-limited, 
linear flow regime. The two are separated by a thermally 
smeared depinning transition^ (not labeled) . (c) Experimen- 
tally obtained domain wall dynamics in layers 1 and 2 in the 
absence of coupling. The two curves cross at H = and 
H = H* . (d) Domain wall dynamics in layer 1 for a coupling 
field, Hi, which reinforces the applied field H, (vi(H + Hi)) 
or works against it (vi(H — Hi)). 

Before considering coupling between domain walls, let 
us consider the simpler case of a single domain wall in 
layer i = 1,2, interacting with a uniformly magnetized 
layer k = 2, 1 (see, for example, Fig. [He)). The interlayer 
coupling J induces an effective coupling field hi, given 
by^ 

hi = mk-rjr- = rrikHi, (2) 

where M s is the saturation magnetization, t is the layer 
thickness, and m = ±1 is the magnetization orientation. 
hi adds to the external field H and also drives the do- 
main wal l 29 i 30 , in turn allowing for a simple experimen- 
tal determination of Hi. To determine Hi, domain wall 
velocities in the hard layer were measured while keeping 
the soft layer magnetically saturated. Through control of 
m 2 and/or H, it was possible to determine wall velocities 
with hi either opposing or reinforcing the applied field. 
We denote these data sets vi(H — Hi) and vi(H + Hi) 
respectively. Plotted in Fig. Hfd), the two data sets are 
separated by 2Hi, allowing a determination of Hi = 120 
Oe and vi(H) [ie. no coupling, see Fig. |2Jc)]. 

The corresponding coupling field and isolated wall dy- 
namics for layer 2 were determined in a different manner. 
H 2 = 220 Oe could be easily found using Eq. which 
gives HiMlti = H 2 M^t 2 (Af s 1,2 are known^ 2 .). Unfor- 
tunately, we were not able to nucleate a domain in the 
soft layer while keeping the hard layer in a single domain 
state and so v 2 (H) had to be measured using a Co(0.5 
nm) layer in a less strongly coupled Pt/Co(0.5 nm)/Pt(4 
nm)/Co(0.8 nm)/Pt filial. 

Now, let's turn to dynamics of coupled walls [Fig.QJb)] . 



The experimental determination of the coupled walls is as 
follows, (i) Two aligned domain walls, at a common po- 
sition xi(0) = x 2 (0), are nucleated, (ii) A magnetic field 
pulse, H, is applied for a time T, under which walls move 
to positions xi(T), x 2 (T). (iii) The new wall positions 
are quasi-statically determined^Si from Kerr microscopy 
images. 

While xi{T) = x 2 (T) for dynamically bound walls, 
xi(T) ^ x 2 (T) for unbound walls since the walls separate 
during their motion. However, the time interval between 
steps (ii) and (iii) is large enough to allow the separated 
walls to relax back to an aligned state under the action 
of effective coupling fields (H = for t > T). Since 
v 2 {H 2 )/vi{Hi) « 10 10 , if xi{T) ^ x 2 {T), pre-imaging 
relaxation of the soft layer wall gives: x 2 ™ s ° d = xi(T). 
Therefore, the experimental technique yields either the 
true bound wall displacement (and subsequently the 
bound velocity) or the hard wall displacement (and there- 
fore the hard layer wall velocity) when the walls are un- 
bound. 

In the unbound state, the hard layer velocity (and 
therefore the experimentally determined velocity of the 
coupled walls), will be that observed for hard layer walls 
under a field H ± Hi [Fig. WIA)] since the walls in the 
two layers are not aligned: +Hi if the hard layer wall 
trails the soft layer wall and —Hi if the hard layer wall 
leads the soft layer wall [Eq. ©]. This is an important 
point, as it allows us to identify the field ranges over 
which v c {H) (the experimentally obtained coupled wall 
velocity) corresponds to unbound motion. The unbound 
(U) and bound (B) regions are labeled in Fig. OHa) in 
which v c (H) is plotted together with vi(H ± Hi) to al- 
low a direct comparison. This allows us to easily locate 
the three critical fields, -ff c i,2,3, which separate bound 
and unbound states [see vertical lines in Figs. [U[a,b)]: 
H c i « 250Oe, H c2 « 750Oe, and H c3 « 1150Oe. 

We first consider the unbound field ranges. In region 
II, H e i < H < H c2 of Fig. da), v 2 (H) > vi(H) [see 
Fig. [UJc)], so that the soft domain wall leads and the 
distance (x 2 — xi) between walls is positive and large. 
The soft wall is so far ahead of the hard wall that the 
latter moves under the action of a positively saturated 
soft layer. When H > H c3 , the situation is reversed: 
vi(H) ^> v 2 (H) [see Fig. J^c) again]. The hard domain 
wall leads and the soft wall is so far behind it that the 
hard wall moves under the action of a negatively satu- 
rated soft layer. A schematic of these regimes is shown 
in Fig. |U 

While we can compare v c (H) to vi(H ± Hi) to ob- 
tain values for the region limits H C j (j = 1,2,3), these 
values can also be evaluated from the experimentally ob- 
tained velocity data in Fig. (He) and the Hi, 2 values. 
Before moving on to analytical and numerical modeling 
results, we explain how this is done using a simple graph- 
ical method. 

In regime II, the walls in each layer move separately 
with the soft wall leading. This can be sustained only if 

v 2 (H-H 2 )>v 1 (H + H 1 ). (3) 
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FIG. 3. (a and b) Experimentally obtained coupled wall ve- 
locity, v c (H) (•) plotted with hard domain velocities in the 
presence of a positively saturated soft layer, vi(H + Hi) (A) 
and a negatively saturated soft layer, vi(H — Hi) (v)- Field 
regions in which bound (B) and unbound (U) coupled wall dy- 
namics are labeled with roman numerals. Vertical solid lines 
represent the region limits (-ff c (i,2,3))- ( c ) v c(H) (•) com- 
pared to hard and soft layer domain wall creep velocities in 
the absence of coupling, vi(H) (♦) and V2(H) (o) respectively. 



Therefore, it is straightforward to define the critical fields 
H c i and H c2 through the equation 



v 2 (H-H 2 ) = v 1 (H + H 1 ). 



(4) 



This equation can be solved graphically using the data in 
Fig. He) to give H°[ w 260Oe and H% « 600Oe, where 
the superscript means that critical field values have been 
determined by the single, isolated, film velocities v\^{H). 
Similarly, in regime IV, the walls are unbound again but 
with the hard wall leading. This can be sustained only if 



Vi(H-Hi) > v 2 (H + H 2 ). 



The equation 



v 2 {H + H 2 ) = vz{H-Hx) 



(5) 



(6) 



now has only one solution, which gives the lower limit of 
regime IV: H^ 3 « 1050Oe. The value H^i compares quite 
well with H c i, obtained from a visual inspection of the 
v c (H) and vi(H±Hi) data above. The bounds H^ 2 , H^ 3 
of region III do not compare so well with H c2 ,H c3 . We 
will comment on that in Sec. IIV CI 

Having considered regions II and IV, we can now turn 
to the remaining regions, regions I and III, which are lo- 
cated around the crossing fields, H = (region I) and 
H = H* (region III). In these field regions, the two walls 



cannot move separately at different speeds, because nei- 
ther Eq. © nor Eq. ([5]) is satisfied. In the following 
Section we argue that in this case a bound state arises, 
for which the common domain wall speed depends on 
Vi(H) in a non-trivial way. 



IV. NUMERICAL AND ANALYTICAL 
RESULTS FOR BOUND STATES 

A. One-dimensional model for wall dynamics and 
numerical results 

In the following we want to introduce a minimal, one- 
dimensional model, which can explain the rising of dy- 
namically bound states and gives quantitative expres- 
sions for the common speed of two coupled walls. Each 
domain wall is approximated by its average position Xi(t), 
i = l,2 [Fig. [Sfa)]. A total field (H + Ht(x)) acts on the 
i— th wall. It is the sum of the external field H and 
the coupling field Hi(x), which depends on the distance 
x — x 2 — x\ between walls. We expect that the coupling 
field Hi is equal to ±ifj, if the two walls are well sepa- 
rated, with the plus (minus) sign applying for the trailing 
(leading) wall. It is useful to make the following general 
assumption for the coupling fields: 



H x {x) = H!f(x) H 2 {x) = -H 2 f(x), 



(7) 



where f(x) is an unspecified odd function, interpolating 
between —1 and +1, as a; varies from negative to positive 
values. Each wall moves with the velocity Vi(H + Hi(x)) . 
A bound state corresponds to motion with 



Vl (H - Hif{x)) = v 2 (H + H 2 f(x)) 



(8) 



for some value x, corresponding to the constant distance 
between walls. If Eq. © has no solution, it means that 
the walls are unbounded (and therefore separated) either 
with the wall in the hard layer leading (v\(H — Hi) > 
v 2 (H + H 2 )) or with the wall in the soft layer leading 
{v 2 (H-H 2 ) >v 1 (H + H 1 )). 

If we define the ratio a — H 2 /H\ between coupling 
fields, we easily find that the solution x — xq of Eq. ©, 

v 1 (H-H 1 f(x ))=v 2 (H + aH 1 f{x )) (9) 

has the form 

H 1 f(x )=G(H,a) (10) 

and the common speed Vb(H) of bound motion is 

v b (H) = vi(H - G(H,a)) = v 2 (H + aG(H,a)). (11) 

Therefore, the specific form of the function f{x) is ir- 
relevant to determine the velocity of bound motion: the 
speed depends only on the external field H and the ratio 
a between coupling fields. Different forms of f(x) give 
different equilibrium distances xq, but the same common 
velocity^. 
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FIG. 4. This figure shows the different field regimes. In the unbound states (thin, black lines), walls in layer 1 and 2 move at 
different speeds, the leading wall at a higher velocity. In the bound states (thick, grey lines), walls move at the same speed. The 
distance between walls is therefore constant. The leading wall is the wall with the higher velocity, in the absence of coupling. 
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pled wall dynamics, vc{H) (•), and the theoretically calcu- 
lated bound state velocity, v b (H) (□). 



where experimental values for and dj are given in Table 
III 





Coupling field 
(Oe) 


exponent 
(Oe 1 / 4 ) 


prefactor 
ln(m/s) 


Hard (1) layer 


Hi = 120 


of = 225.2 


lnu? = 45.3 


Soft (2) layer 


H 2 = 220 


of = 40.1 


\nv% = 10.8 


Bound creep 




ajf « 202 


lnug « hiv'l 



TABLE I. Experimental values for the coupling field magni- 
tudes and parameters for uncoupled domain wall creep dy- 
namics [Eqs. (fl2]l] as well as bound creep [Eq. ([131]. 



We are now going to prove that, in the limit H — > 
0, Eq. (j8|) has a solution which describes a bound creep 
motion such that the common domain wall velocity is 
given by 



v b (H) 



exp 



(13) 



We can now solve Eq. © using experimental data for 
single wall motion (vi(H) and v%{H), Fig. [He)). This 
way, our theory provides the velocity of bound states 
without free parameters. Results are shown in Fig. [S] 
Comparison with experimental data is very satisfying for 
the low field bound state regime, with modest quantita- 
tive agreement in the high field bound state regime. In 
the next Sections we are going to discuss both regimes in 
more detail and derive analytical expressions describing 
the bound dynamics. 



B. The low field bound state regime 

In the creep regime, analytical expressions are available 
for the wall velocities in the uncoupled case [see Eq. (JTJ], 



vi(H) 



v 2 (H) 



vl exp 



exp 



(I)' 



(12a) 



(12b) 



In order to reduce the notation, let us introduce the 
quantity c = f(x), which varies in the interval (—1, +1). 
We have to solve Eq. ©, which using Eqs. (|l"2"j) . can be 
written as 



Vi exp 



H-cHi 



exp 



a 2 



H + cHq 



(14) 

It is clear that walls must move with a positive velocity, 
if the external field H is positive. This requires the sign 
of the total driving fields H ± cHi to be the same as the 
sign of H, which demands that c vanishes in the limit 
H — > 0. Therefore, we use a small H expansion 



c = c H + dH 1 ' 



(15) 



where the value of 7 will be found below, while it is 
straightforward that the leading term is linear. In fact, 
if c vanishes faster than linearly, the coupling would not 
have effect in the limit H — > and a bound state would 
be impossible for small H . On the other hand, if c van- 
ishes slower than linearly, H ± cHi cannot both have the 
same sign as H. In conclusion, using (|15[) we can rewrite 
Eq. (HD as 
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dimensional model and experimental results show only 
modest agreement. Even if our theory correctly an- 
ticipates the existence of a bound state regime around 
H = H* , the agreement between observed (i? C 2,3) an d 
predicted 3) limit field values is not perfect. Fur- 
thermore, the theory (vb(H)) underestimates the exper- 
imental (v c {H)) bound state velocity, Vb(H) < v c (H), 
for H c2 < H < H c3 . Below, we discuss these details 
and, in particular, why the experimental bound state ve- 
locity at H* , v c (H*) ~ 24.5m/s, is significantly larger 
than Ui(iP) = v 2 (H*) = v* ~ 18m/s. In Appendix |B] 
we also give an analytical approximation for the bound 
state velocity in the high field regime. Finally, it is worth 
mentioning that the constant distance xq = X2 — x\ be- 
tween the soft and the hard walls in the bound regime, 
is positive for H < H* and negative for H > H* (see 
Fig. [4]) , because the leading wall in the bound regime is 
the wall with the highest speed in the absence of cou- 
pling. 

Now, let us discuss the disagreement between our the- 
ory and experimental results in the high field bound 
regime. There are three main possibilities to explain this: 
(1) our coupling model is inadequate, (2) the data used 
for v 2 (H) is not representative of the true vi (H) in this 
system or (3) the use of the experimental v\^{H) data is 
not valid for the high field limit. 

(1) In Appendix [XJ we discuss two modifications to 
the coupling: a dipolar coupling (additional g(x) term 
in Eqs. ([7]): see Eqs. (|A1[I ) due to strong stray fields at 
the domain edgea 3 ^ and the use of differing f\(x) and 
f%(x) functions in Eqs. (J7J) (see Eqs. (|A2|I ). However, 
both modifications still lead to Vb(H*) = v* . Further- 
more, since the low field bound regime is well repro- 
duced using only the exchange field, it is questionable 
to make Eqs. ([7]) more complicated. One might also con- 
sider the case in which f(x) is not continuous. For ex- 
ample, we might have a step function, f(x) = — 1 for 
x < and f(x) = +1 for x > 0. This implies that 
the bound state is not characterized by a constant dis- 
tance between walls, but by a continuous interchange be- 
tween the walls. However, this neither solves the issue 
surrounding Vb(H*) ^ v* , nor the discrepancy between 
H c2 . 3 and H* 2 3 . 

(2) As explained earlier, V2(H) was not measured in 
this multilayer but rather in a similar one with an equiv- 
alent Co(0.5 nm) layer. Using this data, we see that in 
the vicinity of H* , walls in layer 2 exhibit flow motion 
wherein V2 = mH with m w 0.022 ms~ 1 Oe _1 . There can 
be some sample to sample variability however and pre- 
vious measurements on a single layer Pt/Co(0.5 nm)/Pt 
film 3 - yielded m w 0.027 ms — 1 Oe _1 . Using this m value 
to model dynamics in layer 2 for the purpose of determin- 
ing Vb{H) at high field changes H* , which is now equal to 
H* = 910Oe. The new m value improves the consistency 
between our calculated H^ 2 and H^ 3 values (700 Oe and 
1070 Oe, respectively) as compared to the experimental 
values, H c2 and H c3 (750 Oe and 1150 Oe, respectively). 
However, the newly calculated value of v% (24.5 m/s) 



remains too low with respect to the experimental value 
v c {H* = 910 Oe) ~ 29m/s 3 4 Note that the film mea- 
sured in Ref. d also had a slightly lower a 1 / 4 value (35.1) 
as compared to a 2 4 (40.1, see Table 1), however this has 
little effect on the predicted bound dynamics since they 
are dominated by the larger a 1 ^ 4 value. 

(3) Finally, our approach, which works well at low field, 
may not actually be appropriate at high field where wall 
dynamics are intrinsically different. At low field, wall 
motion is thermally activated over field-dependent en- 
ergy barriers. In contrast, at high field, wall motion is, 
to a large extent, determined by the internal structure 
of the wall (and associated internal dynamics) 35 ' 36 which 
can actually be modified by interlayer coupling^ 7 -. As 
such, experimentally obtained, isolated single wall veloc- 
ities v\ (H) and v 2 (H) may not be the appropriate build- 
ing blocks to be combined to calculate vj,(H), as we did 
in Eq. (HU). 

V. CONCLUSION 

Exchange coupled Pt/Co layers represent an ideal 
model experimental system in which to study the in- 
teresting problem of coupled interfaces moving through 
physically separate, but coupled, media. Here we have 
detailed the principles behind this system and presented 
both numerical and analytical models of bound domain 
wall motion which compare well with experiment^. Most 
notably, we derive an analytical model with no free pa- 
rameters which describes bound creep. While we have 
concentrated on a one dimensional model we hope out 
results will inspire others to apply micromagnetio 37 ' 38 or 
interface model a 10 ' 39 to this problem. 
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Appendix A: Additional and modified coupling 

A more general expression of Eqs. ([7]) which includes 
dipolar interactions might be 

Hi - H 1 f(x) + D 1 g(x) H 2 = -H 2 f(x)-D 2 g(x) (Al) 

where g(x) = ln[(cP + x 2 )/d 2 ] accounts for dipolar cou- 
pling, d = 3 nm being the separation between the hard 
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and the soft layer. A simple calculation^ can show that 
in the vicinity of the domain walls, dipolar fields can po- 
tentially be larger than H\ 2 . However, as discussed in 
Sec. IIV CI the good agreement between theoretical and 
experimental results at low field suggests that it is H\ 2 
which determine the bound state's stability. 

An alternative generalization of Eqs. is to make 
f(x) different for the two films, 



H 1 = H 1 f 1 (x) H 2 = -H 2 f 2 {x) 
This might mean, e.g., writing 



fi(x) = tanh 



A, 



(A2) 



(A3) 



with Ai ^ A 2 , as is expected for layers with differing 
thicknesses^. 

However, both of these approaches yield Vb(H*) = v* 
since g(0) = and ^(0) = / 2 (0) = 0. ' 



where v* = Vi(H*) and a, = dvi/ <IH\h=h* ■ It can be 
easily shown that the solution x = xq of Eq. (JSj> satisfies 
the relation 



so that 



with 



f(xo) = 



(ai-a 2 )(g-g*) 
a\H\ + a 2 H 2 



v b (H) = v* + a b (H - H*) 



oitt 2 (gi +g 2 ) 
a\Hi + a 2 H 2 



(B2) 



(B3) 



(B4) 



Therefore, in the proximity of H — H* , the common 
speed in the high field bound regime is linear, with a 
slope a b which is in between a\ and a 2 : 



Appendix B: Analytical approximation for the high 
field bound state 

In the high field regime, the walls are no longer in 
the creep regime and Eqs. (Tl2l) cannot be used. Instead, 
we can assume a simple linear approximation^ in the 
proximity of H = H* , 

Vi(H) =v* +(1,(11 -H*), (Bl) 



a 2 < at, < a%. (B5) 

Using the fitting values a\ ~ 0.025 and a 2 ~ 0.12, we 
find a b ~ 0.035, so that 

v b (H) ~ 18 + 0.035(iJ- 850), (B6) 

with H expressed in Oersted and the speed in meters per 
second. 
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